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In biological cells and novel diagnostic devices biochemical receptors need to be sensitive to 
extremely small concentration changes of signaling molecules. The accuracy of such molecular 
signaling is ultimately limited by the counting noise imposed by the thermal diffusion of molecules. 
Many macromolecules and organelles transiently bind to molecular motors and are then actively 
transported. We here show that a random albeit directed delivery of signaling molecules to within a 
typical diffusion distance to the receptor reduces the correlation time of the counting noise, effecting 
an improved sensing precision. The conditions for this active focusing are indeed compatible with 
observations in living cells. Our results are relevant for a better understanding of molecular cellular 
signaling and the design of novel diagnostic devices. 
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Cellular signaling relaying external or internal bio¬ 
chemical cues typically involves low copy numbers of mes¬ 
senger molecules, inevitably effecting appreciable fluctu¬ 
ations in the count of molecular binding events at specific 
receptors A similar limitation by counting noise 

is encountered in modern microscopic diagnostic devices 
to which sensitivity is a key factor [l^. Modern micro¬ 
scopic techniques reveal molecular signaling events and 
underline their inherent stochasticity in living cells [H- 
[T^ . Nevertheless molecular signaling pathways in bio¬ 
logical cells operate at remarkable precision [1^ [l9| . 

The first heuristic argument about noise limitation to 
biological concentration measurements is due to Berg 
and Purcell assuming biochemical receptors to count the 
number N of specific molecules in a volume equal to 
their linear dimension a Q. is then limited by Pois- 
sonian noise SN ~ The time between two in¬ 

dependent measurements is set by the time tjj ^ a^/D 
needed to clear the volume by diffusion, D being the 
molecular diffusivity. Averaging over a time thus 
allows to Am ~ Tm/fD independent measurements, re- 
ducing the noise by the factor Am . The relative accu¬ 
racy to measure a background concentration (c) is thus 
fc^/(c)^ ~ {Da{c)Tm)~^ Hi- When the additional bind¬ 
ing dynamics to the biochemical receptor is explicitly 
taken into account, this relative error becomes [l| 


5c^ 


-b 


(c)2 fcon(c)(l - {n))Tm 7TDa(c)T, 


( 1 ) 


The first contribution stems from the Markovian 
(un)binding to the receptor at detailed balance condi¬ 
tions with binding rate fcon and average receptor occu¬ 
pancy (n). The second term is the diffusional noise, up 
to the factor tt identical to the result by Berg and Pur¬ 
cell [H. The prefactor of the diffusive term in Eq. m 
was recently refined heuristically 11|. Inspired by the 
early ideas of Berg and Purcell a number of studies un¬ 
raveled the crucial role of diffusional noise in biochemi¬ 


cal signaling along with additional features 

such as receptor cooperativity Q and facilitated diffu¬ 
sion [l^. Various experiments suggest that cells indeed 
operate very close to the fundamental accuracy limit [l| . 

Here we extend the result to the case when the 
signaling molecules are not only freely diffusing in the 
cell but actively transported along cellular filaments by 
motor proteins Ml effecting intermittent ballistic ex¬ 
cursions Such an additional active component 

occurs when extracellular signaling molecules are taken 
up by the cell via endocytosis: the molecules are en¬ 
gulfed into submicron lipid vesicle and then intermit¬ 
tently transported through the cell by motors [i^. A 
similar combination of free diffusion and active transport 
occurs when virus particles invade a living cell 27| . How¬ 
ever, even free molecules such as messenger RNA may at¬ 
tach to motors 2^ , or proteins move in directed fashion 
due to cytoplasmic drag [1^. To incorporate the active 
component we employ the theory of random intermittent 
search for hidden targets [s^ which was recently used to 
analyze reaction kinetics in active media (^ . We show 
that active transport enables both faster as well as more 
accurate sensing: an active noise floor exists, but it can 
be significantly lower than the purely diffusive counting 
noise O- This active focusing reduces the noise corre¬ 
lation time and enables the receptor to detect relative 
changes in concentration with higher accuracy. Our re¬ 
sults also have direct implications to the design of active 
components in microscopic synthetic diagnostic systems 
based on molecular signals 11,11 1. 


Model. We consider a signaling particle (vesicle, virus, 
mRNA, or protein) in 3-dimensional cellular space, ran¬ 
domly switching between a passive diffusion phase p with 
diffusivity D and an active ballistic phase a with veloc¬ 
ity v(H) of constant magnitude v = |v| [l^ in the di¬ 
rection of the solid angle H following a Markovian dy¬ 
namics (Fig. [T|). Assuming ideally disordered cytoskele- 
tal filament orientations, the spatial direction of active 
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FIG. 1: Schematic of the model system: a) Each signaling 
particle performs passive thermal diffusion (blue phases) in¬ 
terrupted by active ballistic excursions with constant speed 
and random direction (red phases moving along the black mo¬ 
tor tracks). The duration of both phases is distributed expo¬ 
nentially with mean times Tp^a- When the particle reaches the 
receptor (green sphere) it binds with on-rate kon and dissoci¬ 
ates with rate fcoff. b) Magnification of the receptor region. 


motion events is uniformly distributed. The duration of 
active/passive phases is exponentially distributed with 
mean Ta^p [1^. The concentrations of freely diffusing and 
motor-bound signal particles are Cp(r,t) and Ca(r, I7,t). 
A receptor with radius a is placed at Tq. Then the frac¬ 
tional occupancy n(t) of the receptor by a signal particle 
evolves according to a mean field kinetic scheme obeying 
detailed balance with on/off rates fcon/offi 

dn{t)/dt = fconCp(ro, t)(l - n{t)) - koffn{t). (2) 

Assuming that the particle (un)binds tolfrom) the re¬ 
ceptor only from (to) the passive mode [^, the coupled 
set of equations for the concentrations Ca,p reads 


= DV\ + 

dca{r,n,t) 


—dn - -S{r- ro)——, 

Ta Tp dt 


dt 


= -Vr • (v(f])Ca) - ^ + 


47rT„ 


( 3 ) 


The signaling typically occurs in two stages. In the ini¬ 
tial phase a change in the concentration of the signaling 
particles occurs either by exchange with the extracellular 
space [3 [11 or by variation of the production and/or 
degradation rates. Upon re-equilibration (assumed to be 
much faster than the measurement time Tm) the receptor 
reads out the concentration over the time in the mea¬ 
surement phase. In a diagnostic device equivalent phases 
will be observed after sample immersion and during de¬ 
tection periods. In an optimal signaling setup equilibra¬ 
tion should be as fast as possible while the measurement 
phase should be as precise as possible. We now quantify 
the speed and preeision of the two signaling phases. 


Speed. We assume that the system equilibrates on a 
time scale over which the signaling molecules move a dis¬ 
tance L of the order of the cell size (or that of a cellu¬ 
lar compartment). At this stage we neglect the analyte- 
receptor binding dynamics and adopt a probabilistic view 
of Eqs. dS]). The equilibration time n is then defined by 
the mean squared displacement (MSB), (|r(ri)p) = 

The exact result is (see Supplementary Material, SM) 


(|r^(t)|) = 2 (r;r,)^e--- — 


-I- 3Br„ 


(TT'+rp-i)2 


t(Ta +Tp) 


Tp^{vTa) +^D ^ , ‘iDTp- {vTaf'{l + 2Tp/Ta) 

\ ^ , I' “r 


I+Ta/P 


[l+Tp/TaY 


..(4) 


Over a period of duration Tq -|-Tp, during which the direc¬ 
tional persistence in the active phase causes a nonlinear 
time dependence of (|r^(<)|) and hence a local violation 
of the central limit theorem, an effective diffusive regime 
(|r^(t)|) ~ t is established. Eq. (g]) is a transcendental 
equation for r^, essentially depending on only three pa¬ 
rameters: the typical distance covered in the active and 
passive phases, Xa = vTa and Xp = ^Drp, and the Peclet 
number Pe = Lv/D. To estimate the efficiency of active 
trafficking with respect to diffusion we compare with 
the purely passive equilibration time tq = (6B). The 

results for various Pe values (see Eig. [2}})) typical for 
biological systems are shown in Eig. [2^)-d). 

Active transport is more efficient for larger particles 
(small D) in larger domains, a direct consequence of 
the finite motor velocity and instantaneous direction¬ 
ality of active motion. Namely, in terms of the MSB 
diffusion and active motion display different time scal¬ 
ing (~ t versus ~ t^): considering only pure passive 
and active motion for Pe < 6 diffusion is more effi¬ 
cient. In the intermittent case the motion has a tran¬ 
sient period of duration Ta-\- Ta, which corresponds to a 
parameter-dependent combination of both regimes dur¬ 
ing the relaxation towards the equilibrium partitioning 
between phases a and p. After this transient period an 
effective diffusive regime is established with diffusivity 
Befi = (B -I-a;Q/(3rp))/(l-I-T q/tp), which mayor may not 
be larger than B. can thus be smaller or larger than 
tq. Trafficking of vesicles with B < 10“^^m^/s therefore 
mostly profits from active motion, whereas active motion 
of smaller proteins with B ~ 10/im^/s will only be more 
efficient over large distances as in eukaryotic cells (espe¬ 
cially for neurons), and only if accompanied by signihcant 
phases of passive diffusion. The observed features explain 
why it is profitable for a cell to use active transport for 
trafficking of larger particles BIH: despite deman ding 
more cellular resources. Similarly, active diagnostics [l3| 
can be faster and hence allow for a higher throughput. 

Precision. Since the precision of the receptor mea¬ 
surement of the signal molecule concentration should 
be maximized we consider small deviations from the 
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FIG. 2: a)-d) Equilibration times ro/ri for various Peclet numbers Pe from numerical inversion of Eq. ([4]). Yellow lines denote 
To/Ti = 1. When ro/ri > 1 active trafficking is more efficient. Note the different scales for rolri. e) Pe as function of cell 
(domain) size L and diffusivity D for biological cells [f^ . [3^ . Arrows denote increasing particle size from right (small proteins) 
to left (large vesicles) and cell size from bottom (prokaryotic cells) to top (eukaryotic cells). The border between prokaryotic 
and eukaryotic cells is around fO/rm. For Pe ~ 10"* (vesicles or large proteins in eukaryotic cells), active transport always 
improves the speed. For small proteins in eukaryotic cells and large proteins in prokaryotic cells (10 < Pe < 10^), in addition to 
large Xa a large Xp is necessary for active transport to be more efficient. For Pe < 10“^ active transport is always less efficient. 


equilibrium values n = (n) + Sn and Cp^a = (cp,a) + 
a, taking into account the detailed balance for 
the binding/unbinding transitions via fcon(cp)/fcoff = 
exp{F/[kBT]), where F is the binding free energy (36| . 
Solving Eqs. ([3]) and following the ideas of Ref. Q we use 
the fluctuation-dissipation theorem to obtain the power 
spectrum of the concentration fluctuations. Within the 
linear response regime and for receptor measurement 
times Tm exceeding any correlation times, the relative 
error for the concentration measurement is 

^_ ^{Xa,Xp) 

(Cp)^ ^on(Cp)(l {^')')Xm F(x{Cp^Tpri 

Compared to Ref. 0 for passive diffusion our result dH) 
differs by the dimensionless factor A(a;a, Xp) in the second 
term, for which always A(xo,a;p) < 1 holds (36| . 

We now focus on the transport controlled regime in 
which fcon(cp),/coff > x~^,Tp^ and fcon(cp)/fcoff » 1. 
Hence we neglect the first term in Eq. m and consider 
the remaining active noise floor. If the active excursions 
are short compared to the receptor size, Xa/a <C 1, 

,, , y5aXptanh-i(S;[l + i(W^p)']"") 

(7rxa)2 [i+i(^„/^^)2]i/2 

( 6 ) 


In the biologically more important situation Xa ^ a, 

K{xa,Xp) 1 - [iiXp/a)/{iixp/a), (7) 

which has the lower bound Amin (7rXp/a)^[I—(Trccp/a)^] 
as Xp/a —>■ 0. This might suggest an approach towards an 
infinite absolute precision of the transport term as Xp —>• 
0. However, at fixed total concentration ctot of signaling 
particles in the transport controlled regime we have that 
(cp) = Ctot/(l + Xa/Tp) corresponding to (cp) -i> 0 as 
Xp —?► 0, hence diverging relative fluctuations. Therefore, 
there still exists a noise floor to active sensing but it can 
be significantly reduced as explained below. 

We gauge the sensing precision at equal Ctot and equal 
(cp), corresponding to the lower and upper bounds for 
the gain of active focusing. For equal ctot we compare 
the relative accuracy of measuring different concentra¬ 
tions of freely diffusing molecules in active and passive 
sensing, finding that the precision is always worse for ac¬ 
tive transport as compared to free diffusion (Fig. |3K)-c) 
and asymptotic results in e) and f)) and becomes worse 
with longer active excursions. The reason is that despite 
reducing the absolute fluctuations by active dynamics we 
are measuring smaller and smaller effective concentra¬ 
tions. Conversely, Figs. [3ji)-f) shows that if we compare 
the precision at identical concentration Ctot(l Pxa/Tp) of 
free molecules the accuracy can be improved significantly 
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FIG. 3: a)-c) Fractional variance aa = [(5cp/(cp)^]^^^ of concentration fluctuations at the receptor site for active trafficking 
compared to thermal diffusion, (Tq = [T^DactotTm\~^^'^ as function of the average distance traveled in both phases in units of 
the receptor radius a. In a)-c) the total concentration ctot is kept constant, hence decreased sensing accuracy is solely due to 
exceedingly small equilibrium concentrations in the passive phase, d) CTa/no for equal equilibrium concentrations (sc in parts e), 
f)) in the passive phase (cp) for different total concentrations = CtotCl + Ta/rp): a strong increase in sensing accuracy is observed 
for Xp/a —>■ 0. a)-d) use the full results ([5]) and (IS7II in SM [^ . e) Horizontal cross-sections of (Ta/cro at large and small Xp 
(symbols) compared to the approximations (|6]l, 0. f) Vertical cross-sections at large and small Xa (symbols) compared to the 
approximations ®, 0. Note the different scales for (Ja/(Jo and the excellent agreement between full and approximate results. 
The non-monotonicity at small Xp is due to the interplay of signal focusing (smaller A) and decreasing (cp). 


as long as Xp < a (otherwise, in this regime active traf¬ 
ficking does not affect sensing precision). The ingredi¬ 
ents necessary to understand this reduced counting noise 
are: (i) only fluctuations on a scale ~ a are relevant for 
the sensing accuracy and (ii) the receptor only ’sees’ free 
particles. Hence, at finite temperatures perfect signaling 
corresponds to the situation in which upon release from 
the receptor the particles immediately bind to a motor 
and are swept away over distances > a. Concurrently 
diffusive displacements must be < a to assure focused 
delivery in the sense that any unbinding form the motor 
only contributes if it occurs at |r — ro| ;< a. This re¬ 
duction of local concentration fluctuations by means of 
intermittent active excursions is exactly the active signal 
focusing mentioned above. 

Conclusion. The last years have seen significant ac¬ 
tivity to egjlore the counting noise for purely diffusive 
scenarios [iMllj and to analyze the speedup of recep¬ 
tor binding due to the topological coupling of one- and 
three-dimensional diffusion for gene regulation in the fa¬ 
cilitated diffusion model [3 El. In this Letter we fill the 
apparent gap in the quantitative assessment of the count¬ 
ing noise experienced by biochemical receptors measur¬ 
ing the local concentration of signaling molecules or com¬ 


pounds in the case when an active transport component 
is present. This occurs for various signaling particles 
(proteins, mRNA molecules, vesicles containing signaling 
molecules, or viruses) by direct shuttling of these signal¬ 
ing cues by molecular motors or by cytoplasmic drag. 


Compared to the purely diffusive signaling considered 
so far we showed that the counting noise (the limit to 
accurate receptor measurement) for active sensing can 
become significantly reduced due to active focusing. The 
only contributions to the counting noise stem from par¬ 
ticles, which are actively transported to within the parti¬ 
cle’s typical free diffusion distance to the receptor. This 
reduces the correlation time of the receptor occupancy 
noise and renders the averaging over a measurement time 
Tm more efficient. The importance of active signaling in 
cellular regulation is well recognized (3 |. In agreement 
with our results, in biological systems active transport is 
indeed employed to move larger particles (e.g., vesicles or 
viruses) with intrinsically small D 3 3-24, 27, 35|. As 
a result even during longer periods of detachment from 
motors these particles barely move 3 3- 2^ 35 L The 


typical experimental values Xq ~ 0.5 — 10/rm [2.31-125 


for a ~ 1 — lOnin in fact fulfill the requirements of our 
model for signal focusing. However, as discussed here also 
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smaller particles such as mRNA and proteins experience 
active motion components [H, , effecting active focus¬ 

ing for their detection. In living cells the motor tracks are 
often not ideally disordered, as assumed here, but biased 
towards the receptor [S^. An expected net directional 


bias towards the receptor, while not imp eding signal fo¬ 
cusing as long as Xa ^ a and Xp < a (33| , would enhance 
the rate of delivery and simultaneously increase the local 
concentration (cp(ro)) at the receptor for equal Ctot- Sig¬ 
nal focusing is thus inherent to active cellular signaling. 
Conversely, despite the great technological advance over 
the past years, molecular motor-powered diagnostic de¬ 
vices have not yet demonstrated a performance beyond 
the existing passive techniques, but a large superiority is 
much anticipated 0. Our results confirm these expec¬ 
tations and present a first rigorous theoretical basis for 
their systematic improvement and development. 
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Supplementary material 
Signal focusing through active transport 

Aljaz Godec^’^ and Ralf Metzler^’^ 

^Institute of Physics & Astronomy, University of Potsdam, D-lff76 Potsdam-Golm, Germany 
^National Institute of Ghemistry, 1000 Ljubljana, Slovenia 
^Department of Physics, Tampere University of Technology, FI-33101 Tampere, Finland 

In this Supplemental Material we summarise the calculations and full results for active transport coupled to reversible 
binding to a receptor. 


SPEED 


In this part we neglect the binding to the receptor (see main text for argumentation) and invoke a probabilistic 
interpretation of the model (i.e., concentrations —^ probability densities). By Laplace transforming Eqs. (2) and (3) 
in the main text in time and Fourier transforming in space, (j,t) —> (k, s) can be solved for Ca and Cp. The Laplace 
image of the mean squared displacement (MSD) for a particle starting at the origin in the passive phase is obtained 
from (|r^(s)|) = —V^(/ Ca(in + Cp)|k=o and interchanging the order of differentiation with respect to k and integration 
with respect to LI: 


(Ir^WI) 


6D 




+ Tp- (TaTp)-V(S + Ta ) 


(SI) 


which can be inverted exactly to give Eq. (4) in the main text. The equilibration time Ti is defined implicitly by 
(|r^(ri)|) = L^, L being the linear dimension of the cell or its compartment. The transcendental equation for Ti is 
solved numerically. 


SENSING PRECISION 


The Fourier transformed linearized Eqs. (1) read 

—iujSn{uj) = —T~^Sh{io)kon{l — {n))6cp{ro,io)koe{n)f)6F{uj) 
—ioj5cai}f;Ll,uj) = —(iv(n) • k + T”^)5c~a(k, 17, w) + (47rTp)“^(5cp(k, w) 


—iuiScpCk, uj) 


-{Dk^ + Tp i)(5cp(k,w) + r„ 1 


Sca(k,I},uj)dLliujSfi{uj)e 


where we have used the constraint imposed by detailed balance Skon/kon — dkos/kos = SF/^ksT) when linearizing 
the first of Eqs. (1) in the main text. Solving for dn{uj) we obtain 


5h{uj) 


_ kos{n)l36F{uj) _ 

Tc“^ - zw (l + fcon[l - («)](27r)-3 / dkS(k, w)) ’ 


(S2) 


where j3 = l/iksT) and we have defined the correlation time = {kon{cp) — kos) ^ for two-state Markovian switching 
[l[ as well as the two auxiliary functions defined as 

S(k,a;) = -ioj -b Dp -b - (47rTpTa)“^'I'(k, w) (S3) 


'I'(k, w) 


/■2Tr 


dO sin(0) 


dip 

i{uj — v(0, p) ■ k) 


(S4) 


We are interested in the low-frequency limit and hence take 'I'(k, w) ~ 'I'(k, 0). In this limit the integral over p in 
'I'(k, w) can be evaluated as a contour integral along the unit circle via the method of residues, while the second one 
is amenable directly leading to 


47r 

4^(k, 0) = — tan 
vk 


-1 


VTak) , 


(S5) 
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where k = |k|. Using Eq. (IS5I) we can obtain 5h(ui) as a function of 5F[uj) explicitly. The linear response function of 
the receptor occupancy (the coordinate) to a change in the free energy difference between bound and freely diffusing 
species (the thermodynamically conjugate force) is then Sn{uj)/SF{u!) and is related to the power spectrum of n via 
the fluctuation-dissipation theorem 5'„(w) = 2/(/3w)Im[5h(u;)/(5F(w)]. A change in concentration is equivalent to a 
change in F as a consequence to a change in chemical potential, 6cp/{cp) = j35F. Using this it can be shown [i| that 
S'c(w) = {l3{cp)Y\5n{bj)/5F(u:)\~‘^Sn(!^) leading to 


>S'c(w) 


ojksT Sn{uj) 


(S6) 


We assume that the receptor measures and averages the concentration over a time much longer than any correlation 
time. We are thus interested in the low frequency limit Sc{oj —>■ 0) leading to 5cp = Sc{uj —>■ 0)/rm and hence take 
S(k, w) Ri S(k, 0). Finally, cutting the integral in the inverse Fourier transform |k|j„aa; = 7r/a to correct for a finite 
target size (and thereby avoid the UV divergence) we arrive at Eq. (5) in the main text with 

HXa,Xp)= f 

Jo V \TTXpqJ 

The integral in Eq. (IS7I) can be evaluated exactly in the limits when the active excursions are either very short 
Xa/a <C 1 or very long Xa/a ^ 1 compared to the size of the target and gives Eqs. (6) and (7) in the main text. 

Now we have to evaluate the equilibrium concentration in the passive phase in the general case, where signaling is 
not necessarily transport controlled. In absence of a directional bias the concentration will be uniform in both the 
active and passive phase. Hence we only need to solve a system of linear equations for (n), (cp) and f (ca)dfl. It can 
be shown that the steady state concentration of freely diffusing molecules (cp) is given by 


1 - 


tan (TrXaq/a) 


^n/r 


dq. 


(S7) 


/ \ _ ^tot ^ ^off 

2(1-|-Ta/Tp) 2fcon 

where V is the volume of the cell or its domain. For fcofi/(fcon(cp)) —)> 0 we obtain (cp) = Ctot/(l + TajTp), where for 
convenience we have absorbed the term 1/U in the total concentration ctot — V~^ —>■ Ctot- Expressing the ratio Ta/xp 
in terms of Xp and Xa, Ta/xp = we observe that (cp) and hence the relative accuracy along with Xa and Xp in 

fact depends on the factor as well. 


/ Ctot - V-^ 

^off 

\ 1 ^tot^off 

V2(l -1- Xa/Xp) 

2^011 / 

/ kon{i+Xa/Xp) 




(S8) 
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